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An asymptotic expression is obtained for the current distribution on the outside sur- 
face of an infinitely long, perfectly conducting, hollow cylindrical antenna that is fed by an 
infinitesimally narrow circumferential gap. This asymptotic expression involves two 
series. The first series is expressed in reciprocal powers of log (2\z\ljr 2 ka 2 ), where \z\ is 
the distance from the gap, log r is Euler's constant, k is the propagation constant, and a 
is the radius of the antenna. The second series is a similar series multiplied by l/(k\z\). 
The first series is dominant and its first five terms yield values for the magnitude and phase 
of the current that for even moderately thick antennas (circumferences as large as X/3) 
are accurate to within about one percent in as close as X/3 of the gap. This is shown by 
a comparison of the values of the current obtained from these terms with the numerically 
computed values of Duncan [1962]. Asymptotic expressions for the current found in the 
literature resemble the first term of this dominant series and are accurate only at rela- 
tively large distances from the gap) — except for very thin antennas. 

1. Introduction 

Consider the usual model of an infinite antenna composed of a hollow circular cylinder 
with its axis lying along the £-axis of a cylindrical coordinate system with coordinates (p, 0, z). 
The cylinder is cut in two at 2=0 and the antenna is assumed to be excited by a sinusoidal 
generator producing an oscillating electric field E g e iut (that is independent of </>) across the 
gap in the cylinder formed by the cut. The voltage across the gap is 

V=- f Edz 



--/. 



which may be held constant while the gap width is allowed to approach zero. As is well 
known this limiting process of allowing the gap width to go to zero leads to a singularity of 
the quadrature component of the current I(z) on the outer surface of the cylinder. The 
effect of such a gap has been considered, among others, by Infeld [1947], King [1956], Wu 
and King [1959], Chen and Keller [1962], and Duncan [1962]. 

The current on the outer surface of the cylinder is given by the formula * 



7( * )= id- 



where 

Jc=2t\- 1 

X— wavelength in free space 

d=radius of antenna 

Zo=-yfiJJe 
€= dielectric constant of free space 
H= permeability of free space 
0=(F-7 2 )* 
H^Q3a)=J Q (fia)-jT (fia) 

> See formula 13 of Duncan [1962]. 
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H?($a)e- m *dy . 
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BRANCH CUT No. 2 




Figure 1. Contours of integration in the y-plane. 



BRANCH CUT No. I 

y-PLANE 



The last two functions are Hankel functions of the second kind, the time factor having been 
taken as e jcct . The path of integration in (1) is along C in figure 1 and is thus along the real 
axis except at y= — k where it is downwardly indented and at y=k where it is upwardly indented 
to avoid the branch points of fi at y = ± k. The angle of /3, written Z/3, is chosen so that Z/3=0 
for —k<^y<CJc. 

Duncan [1962] has recently derived expressions for I(z) that permit its computation for 
very small values of z and has obtained accurate plots of the real and imaginary parts of I(z) 
as one approaches the singularity at 2=0. While this approach also makes possible accurate 
values of I(z) out to f =k\z\ of 5 or 10, the work involved increases and the accuracy decreases 
as f becomes large. 

This paper develops an asymptotic series for the current 



I(z)5±-y-e 

where A n and B n are constants and 



■*r 



lift A„ j<& Bn "I 

\J={ (logpf)" f ^2 (l0gK)"J 



'jr 2 k 2 a 2 



(2) 



(3) 



It is shown that the first five terms of the first series yield a value of I(z) that is remarkably 
accurate even for rather small values of $=k\z\ and for thick antennas. Thus, for example 
(see figs. 2 and 3), \I(z)\ is given within an accuracy of 1 percent for all antennas thinner than 
a = 27ra\- 1 =0.30 for all f=27r|s|\- 1 > 1.8 (approx.). 

2. Modification of the Integral for I(z) by Contour Integration 

This integral for I{z) given in (1) may, by reference to figure 1, be written 



7(g) = 



kaV 



Lira 



!' 



jZ € ^ Jc 



pHPtfia) 



■jy\z\ 



-dy. 



(4) 



By distorting the contour C between A and B into C, one sees that C in (4) can be replaced 
by C . It can readily be shown that in the limit of e— >0. and R— >oo only the portions of the path 
C lying along branch cut 1 contribute to the integral, and hence 



jkav r*-i-rm 



'(/3a) 



H^iPae- 1 *) 



PH^ (jSa) 0e ~ *Hf> (fiae - ir ) 



} 



— iY'zl 



dy. 



(5) 



Decreasing the angle of the argument fia of the Hankel functions by tt leads to the following 
equations : 



418 



HP ipae-f') = -Jo(0a) -jY t tfa) = -flj 1 ' G3o) 
fl?> Ctae"*') = J, (rfto) +iF, (/3a) =ff«> (0a) . 

These together with the substitution y=k(l—jy) and the well-known identity 

2 



j 1 (x)r (^)-j (^)r 1 (^)= 



■KX 



makes it possible to write (5) in the form 

jWe 
where 



T(~\ J We " H C e — 

1{Z) 7rZ Jo „VS0»«)- 



+*%«)] 



£=k\z\, a=ka, and v=^y(y+^j)- 



(6) 

(7) 

(8) 
(9) 



The change in the contour of integration thus leaves us with a more manageable integral for 
large values of f in that only the small argument behavior of the Bessel functions is important. 

3. More Tractable Expression for the Asymptotic Behavior of I(z) 

We may write (8) in form 



mJv^T 



where 
and 



7rZ 

T 2j p e-*du 

f Jo V 2 [JIM+Yl 



(v<*)\ 



r)=jju(u+2ft). 



(10) 

(ID 
(12) 



Within the region Irja^b, where b is a small positive number, the Bessel functions may be 
replaced by their small argument approximations. This amounts to setting 



l + €i(r?a) 



dm J 2 (va)+Yi( V a) 1+ /2 log r^y 



(13) 



and neglecting e t (rja). This leads to a maximum relative error in the approximation that is a 
function of b. For 6=0.5 this maximum error is about 1 percent. 
Letting Ui be the value of u for which 



we have 
where 



\rj a \==-^u 1 \u 1 + 2j^\=b 



(14) 
(15) 
(16) 



It is therefore seen that for any choice of a and b one may make U\ as large as one chooses by 
taking f large enough. A tabulation of Cfor 6=0.5 and various values of a is given in table 1. 



Table 1 



a. 


C for 6=0.5 


T 


t-1 


0.066r-i 


Ul 


0.02 


24.96 


6.34X10-*^ 


1576^ 


104f 


2r 


0.08 


6.09 


1.02X10-2^-1 


98. 5r 


6.50r 


2r 


0.30 


1.19 


0. 143$^i 


7.oor 


0. 462f 


0. 462r 


0.60__ 


0.34 


.571J-1 


1.75f 


.lier 


.ner 


1.20 


.09 


2. 28J-1 


0.44^ 


.029f 


.029r 
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From (11) and (13) 

f Jo ,,r 1+ ^2 log r5«\n T f J., ,*z?(i,«) 

Since it is possible to choose an A such that \D(rja)\^> (Aalr)])' 1 for all |i?a|>&, we may make 
use of (12) to obtain the inequality 



Ur^q^ r€i du<Aa TY (18) 



which, using (15), can be written 



If 



£s»fv&~ <«> 



Thus with the neglect of e t in (17) and exponential term in f we may write 



" 2 [ 1+ 1O: 



g ^) 2 ] 



(20) 



where M is % or any quantity very large compared to 1 . Since the integrand of the first integral 
in (17) has nearly a constant angle, the percentage error in dropping the integral involving 
€1 is certainly not much more than the maximum error in dropping e x in (13). 

Expansion in (20) of [(u/£)+2j]~ 1 by the binomial theorem leads to a series in reciprocal 
powers of f 



where 

"~ Jo i+(>g^- 2 ) 2 ' (22) 

At this point to insure uniform convergence of the binomial series, M is to be taken slightly 
smaller than u 2 , where u 2 is the smaller of U\ and 2f. 
From (12) 

where 



V 2 =v 2 o~v$ (23) 



^Vr (24) 

Consider now a function/ of rj 2 in the range in which w<<f. By a Taylor expansion 

f(v 2 )^m)-\ri^- (25) 



Since, however, 
(25) may be written 



dr\l drjldu 2 du 
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jv? df(r,l) ^ 



Applying (26) to (22), we have 

The application of (26) in the outer part of the range <u<M will not be very accurate, but 
since the integrand drops off as e" u , one may ignore this. Integrating the last integral by parts 
and neglecting the term containing e~ M in the integrated part, we have 

T n ^Tl+^[(n+l)T° n+1 -Tl +2 ] (27) 

where the superscript zero indicates that t] is to be replaced by t? . 

If in (21) one neglects terms of (l|f) 2 and higher, then one may write that 



By (27) this becomes 



Tq±TI+^-[2TX-T%\. (28) 



Equations (10) and (28) provide us with an asymptotic expression for I(z) involving 
integrals of the form 



J° l+(±log M r/) 2 



(29) 
where use has been made of (24) and where 



2f 

4. Asymptotic Series for I{z) 
Integrating (29) by parts using the function 



(30) 



7T 

^4(x)=^+arctan x (31) 



we have 



i T Q n =S u n e- u A (hog ur j)T+ J {n n -nu n - l )e~ u A ^log urj) du. 

Since Ay- logurj )->0 as u — >0 and Mis large, this reduces to 

T°»=W H -nW n -i+M n e- M A (hogMrj) 

^Wn-nWn-t (32) 

W w =7r| u n e- u A(-logurj\du. (33) 



where 
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Collecting our results thus far, we have from (10), (28), and (32) that 

I{z)J^~^ \Wo+^ (4W 1 -2W a -W 2 )^ (34) 

Corresponding to the expansion of the arctangent in reciprocal powers of the argument, 
one has 

/l \ (— l) fc+1 /l \-(2*+D 

A {^ UT 0=Smr\^ ogUTJ ) • (35) 

This series converges for 

1 

>i 



and thus will converge when 



-log Urj 



-log^r> : y7ror^<0.066T- 1 . (36) 



Let Uz be the smallest of the quantities w 1 = (7f, 2f, and 0.066r l then (35) converges 
uniformly for 0<u<M<Cuz, and with neglect of terms of order e~ M 



with 



lnk ~ Jo (log^j) 2 ^ ^ 



This follows from the fact that a uniformly convergent series may be integrated term by term. 
It is here understood that M is to be taken only negligibly smaller than u 3 . It is seen from 
table 1 that u 3 , and hence M, can be taken as large as one desires by increasing f. 
Equation (38) may be written 

— 1 C M u n e~ u du 

^"(logTJ^Jo r i+ log^T* +1 ' (39) 

L l0g TJJ 

Using Taylor's series with a remainder term to expand the denominator, we have (see Widder 
[1961]) 



where {s} as a superscript is a pseudoexponent denned by the equation 

r (2£+l)(2&+2) . . . (2fc+«), s>0 

1, s=0 (41) 

J(2jfc)(2i-1) . . . (2k+s + l)]-\ s<0 



(2k+iy s i--= + 

and 

Letting t— —log v/log rj in (42), one may write 



-lo^r -iog^ T 

logrj L lOg TJ J 

* N\ Jo (l-f)2*+^+2 



. J2^+l)^+i) Cl^j Llogrj 'J (42) 

UN ATI I /-i j\ 9fc_l- A7J-9 ^6. 



viV 



_ (2fc+l)'* +1> f " ( log J . 

* AH (log TJ) y+1 J t / log yT jY +JY+2 p" ^ ; 

\logrj7 
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On substituting from (40) into (39), one obtains the equation 

/ _ -* rf (2fc+o"» q« +B > (w) i (44) 

'- (logrj)^ 1 ^ S\ (l0gri)' +jKwW J (44j 

where 

/•Af 

C^=(-1) S u n 6- M (log^)V/u (45) 

and 

R' N (n)= u n e~ u R N du. (46) 

Jo 

It will be shown later in section 5 that as £— >oo and hence t-^0 the remainder term R f N (n) 
approaches zero as |log r7|~ (2V+1) . This insures that (44) is a legitimate asymptotic expansion 

Of Ink- 
It is convenient to replace (45) by the equation 

C ns ={— l) s \ u n e- u (\ogu) s du (47) 

which involves again only the neglect of a term proportional to e~ M and to write (44) in the form 

-1 jg (2fc+l)M On. m 

lnk = <\0g rjr+ l U S\ (logrj) s ^ 

to indicate an asymptotic series that must be terminated at some finite but unspecified term. 
The remainder term is then very roughly indicated by the first neglected term. A precise 
determination of the overall error is made later in section 7 by comparing the results obtained 
from our asymptotic series with those obtained numerically by Duncan [1962]. 
Integrating (47) by parts, we have 

(-iyC n -Y u n+i e -u (logj/) s+l y + j_ [™ [ U "+ l -(n+\)u n ]e-\\oguY+ l du 

i__ § i 1 Jo 8~t 1 J o 

or 

Cn+l, s+l~(^~l~l)£x s+1 (s-f"l)Cw f «• 

Thus 

On. s =nC n - 1> s—sCn-!, a - t (49) 

and, in particular, 

C/i s =Co s 5Co i5 -i Ww 

=2C„ s -3sC , s - 1 +s(s-l)(7o, s - 2 . (51) 

From (37) and (48), since in (37) W n can be approximated by a truncated series, 

Wo+i Wx-2Wo-W a ] 
2j 

" X S 2H-1 VlogriJ § s! (logri) s L a+ 2f (4 ^ - 2C °- C *>>] (° 2 > 

and thus, finally, by (34), (41), (50), and (51) 

2irVfi-^ (a) ( ") C— l)* +1 C2yt4-2') !s - 1 )7r 2 *r i ~1 

/(2)S ^S S ( idog rfi™ [ C >-i (^.-t+*(-l)^)] (53) 

<7 s = C 0s =(-l)* f V*(log u) 8 rfu. (54) 
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where 



Equation (53) may be written 

2wVe -' risi 



K*h 



■[-(<») 
L»=i 






B n 



A 



Z L^l (log K ) " ' f £3 (log pf )- 
where -4's and B's are constants and —log rj is replaced by log p£ ; thus 

-2i 



p= 



r 2 « 2 



(55) 



(56) 



These are the same equations as (2) and (3) in the introduction. Vainshtein [1959] has obtained 
a similar type of expansion for a related integral. 

The second series in (55) is not, strictly speaking, a part of an asymptotic expansion of 
I(z) since terms in 1/f fall off faster than any power of (log pf) _1 . It can, however, be used 
to increase the range of f over which (55) is a valid approximation. 

By direct comparison of (53) and (55) we see that 



A\ — 6o 




^ 2 =-a 


B,=\ C 


7T 2 
^3= — 77 60+62 


B 3 =-(G +C 1 ) 


At^a-C* 


2 O 

Z 4 = — — 6o+36i+-- 62 



(57) 



^5=t C -2x*C 2 +C t B b =2Tr\C +C 1 )-QC 2 -2C 3 
5 

5. Magnitude of the Remainder Term R' N {n) 
If u is in the range 0<u<r"" 1/2 , then from (43) 



i»»i< 



(2fe+l) |y+11 
■JVl|logrir+ I 



logT+|i 



1, , x . 

2 l0gr +2-? 



2/L-+AT+2 



L N+i J, 



< 



(2fc+l)W-ii2 2 *+"+ 2 |logul JV + 1 



(iV+l)!|logrir+ 1 



(58) 



Since |log Vrj\ in (43) for real v is always greater than t/2, we have 

^ (2^+l) |N,+li |logrj| 2,t+1 llogwr +1 

| ft AT | <. / \ 2*+2V+ 2 



(tf+l)l(p 



(59) 
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even outside the above range. 

These inequalities applied to (46) yield an upper bound for R' N (n) expressed by 

+ l l^^' f ^-"i lQ g <+'<*«}• (60) 

w 

The second integral approaches zero as e~ T ~" 2 and thus for small r (large f ) is negligible. One 
sees therefore that B' N (n) goes to zero as 

constant 



|log rj\ N + l 
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(61) 



6. Computation of the C n 

From (54) 

<? w =(-l) n J o " e-"(log u) n du=P n +(-l) n Q n (62) 

where 

P„=(— 1)» fV*(log^)\/<>0 (63) 

J o 
and 

Q„= J e~ u (log u) n du>d- (64) 

In (63) let x= — log w; then 



P n = e- e ™<- x) x n e- x dx 



^o p! Jo ^op!(p+l) B+l 

and one obtains a convenient alternating series for P n . 
To obtain the Q n we set 

&= ^ M ^"(logu) w ^+e w (M) (60) 

where the truncation error e n {M) in using M rather than infinity as the upper limit is given by 

^(M)=J^- w (logi/)\/^ (67) 

(x \ X 
l+Tf)<- A f for all x>0 and hence e n (M) is less than 

E tl {M) = j^^\-\M\ogM+xyrx 
which on successive integration by parts reduces to 
E (M)= e —m± {mo * M)k 

= ao g M)-e-[l +M ^+^=^+ • • • + IM] ^ W ] (68) 

Setting n=o and]M=16, we have 

J B 5 (16) = (2.77) 5 6- 16 [1.124]-=2.0SX10- 5 . 
For smaller n the error is even less. Thus, we may with negligible error set 

Qn= f e- u (\ogu) n du, n=0, 1,2, . . . (69) 

and evaluate these integrals by numerical quadrature. Actually, one knows that 

C Q = C e~ u du=l (70) 

Jo 

and from standard integral tables that 

0i= e~ u log udu=y = Eulev's constant. (71) 

Hence numerical integration is not required in (69) for n=0 and n=\. 
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A table of the C n and the A n and B n obtained from them by (57) is given below: 

Table 2 



n 


C n 


A n 


B n 




1 

2 

3... 


1.0000 
0. 5772 
1.9781 

5. 4449 
23. 561 
117.84 


1. 0000 

-0. 5772 

-1.3118 

0. 2520 

3. 9969 


0. 5000 
-1.5772 


4 

5 


-0. 2360 
8. 3743 



7. Accuracy of Results Obtained Using Asymptotic Series 

Since several approximations were made in deriving; the asymptotic series (55) for the 
current on an infinite cylindrical antenna, it is important to check the results against the very 
accurate results obtained by Duncan [1962] by numerical integration. For this purpose only 
the first five terms of the A coefficient series were used. The computations were made for rather 
thick antennas, a=0.02, 0.08, 0.30, 0.60, and 1.20, and some of the results are shown in 
figures 2 to 10. 
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Figure 2. Magnitude of the current for a = 0.0S. 
"Asymptotic" refers to first five terms of the Aperies in (55). 
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Figure 3. Magnitude of the current for a=0.30. 
"Asymptotic" refers to first five terms of the A series in (.55). 
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Figure 4. Magnitude of the current for a=0.60. 
"Asymptotic" refers to first five terms of the A series'in (-55). 
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Figure 5. Magnitude of the current for a= 1.20. 
"Asymptotic" refers to first five terms of the A series in (5.5). 



The results for a =0.02 are not given since at a =0.08 (see figs. 2,6, and 9) the results are 
already quite accurate even as close as a quarter wavelength from the gap. For thinner 
antennas the agreement would be even better. One observes (see figs. 4, 7, and 10) that for 
o:=0.60 the asymptotic series begins to give accurate results only when f is greater than 3.5 
or 4. For a circumference approaching one wavelength (a=l) the asymptotic expansion 
ceases to be even approximately valid for determining the magnitude of the current as can be 
seen from figure 5 (a=1.2). Surprisingly enough figure 8 shows that the phase angle of the 
current for o=1.2 is somewhat better predicted, being less than 20 deg off from the true value. 
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680071—63- 
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Figure 6. Phase angle of the current for a=0.08. 
"Asymptotic" refers to first five terms of the A series in (5.5). 
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Figure 7. Phase angle of the current for a=0.60. 
"Asymptotic" refers to first five terms of the A series in (55). 
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Figure 8. Phase angle of the current for <x=1.20. 
"Asymptotic" refers to first five terms of the A series in (55). 




Figure 9. Imaginary part of the current for a=0.08. 
"Asymptotic" refers to first five terms of the A series in (55). 
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Figure 10. Real part of the current for a=0.60. 
"Asymptotic" refers to first five terms of the A series in (55). 



The first term of the asymptotic expansion (55) taken alone yields the asymptotic formula 



J( 2 )S 



2tV 



-jk\z\ 



Z « log -^ 



(72) 



T 2 ka 2 



which is accurate for antennas of the thickness here considered only for rather large £=k\z\. 
This is shown in figures 9 and 10. Even for the a =0.08 case the plot of the first term fails to 
give close agreement even at a f of the order of 10. These two figures illustrate, in addition, 
the degree to which the first five terms give the real and imaginary parts of the current. Since 
the accuracy is approximately the same for each part, we have given only the imaginary part 
for a =0.08 and only the real part for a =0.60. In the latter, of course, the agreement is not 
too good. 

Formulas similar to (72) have been obtained by others. Chen and Keller [1962] obtain the 
same formula but with the argument of the log the negative of that in (72). Northover [1958] 
derives a formula in which the argument of the log is z/2ka 2 . Hallen's [1956] formula agrees 
with (72) except that r is to the first power. All of these formulas are correct in that they lead 
to the right functional behavior for f — >oo . It has been checked, on the other hand, that in 
the case a =0.08, at least, formula (72) is more accurate for finite f than any of these variations 
of it. 

We have also determined that for a=0.08 and 0.60 the error remaining after taking the 
first five terms of the A coefficient series is given for intermediate f very nearly by the first 
two or three terms of the B coefficient series of (55). 
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